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^ ■ Abstract. We carried out numerical calculations of a contribution of the waterfall field to the 

primordial curvature perturbation (on uniform density hypersurfaces) (, which is produced 



during waterfall transition in hybrid inflation scenario. The calculation is performed for a 
^ ■ broad interval of values of the model parameters. We show that there is a strong growth of 

'nI" , amplitudes of the curvature perturbation spectrum in the limit when the bare mass-squared 

of the waterfall field becomes comparable with the square of Hubble parameter. We show 
^f) • that in this limit the primordial black hole constraints on the curvature perturbations must 

r~^ . be taken into account. It is shown that, in the same limit, peak values of the curvature 

^D I perturbation spectra are far beyond horizon, and the spectra are strongly non-Gaussian. 
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1 Introduction 

In last two years an interest in hybrid inflation models was very much revived [1-8] (see also 
earlier works [9, 10]). The main questions which were discussed are dynamics of the waterfall 
field and its influence on the total spectrum of density perturbations produced by inflationary 
expansion. It was suggested, in particular, that fluctuations in the waterfall field could lead 
to non-Gaussian curvature perturbation at rather large scales (even, possibly, at cosmological 
scales). Note that, in general, hybrid inflation models remain theoretically attractive up to 
now, especially in the context of supergravity and string theories (see, e.g., [11, 12]). 

In many cosmological scenarios the period just after the end of inflation, i.e., the inflaton 
decay and the subsequent evolution of the decay products to a thermal equilibrium starts 
with preheating [13-16]. One of the most studied models is hybrid inflation with tachyonic 
preheating. It had been shown in [17] that preheating after hybrid inflation goes through 
the tachyonic amplification due to the dynamical symmetry breaking, when one of the fields 
rolls to the minimum through the region where its effective mass-squared is negative. In the 
process of this rolling amplitudes of field fluctuations grow exponentially leading to a fast 
decrease of a height of the inflationary potential ( "false vacuum decay" ) . 

As is shown in the recent papers [1-4, 6-8] the power spectrum of the comoving curvature 
perturbations from the waterfall field in hybrid inflation with tachyonic preheating is very 
blue: it depends on comoving wavenumber k like (k/k^)"^ (for k < k^) and is negligible on 
the cosmological scales. Absolute value of the spectrum amplitude at /c ~ /c,, (as well as the 
value of /c*) depends on model parameters. In principle, this spectrum at fc ~ /c* may be 
quite substantial, and, in this case, it may be constrained by data of primordial black hole 
(PBH) and relict gravitational wave (GW) searches. 

The amplification of curvature perturbations after preheating (e.g., in two-field infla- 
tion models) and blue spectra of the curvature perturbations on the constant energy density 
hypersurfaces, C; behaving like ~ /c^ on super-horizon scales (in a case of the quadratic infla- 
tionary potential) had been predicted in many works (see, e.g., [18-20]). Another example 
where the steeply blue (~ k^) curvature perturbation spectrum is predicted is one of mod- 
els of false vacuum inflation, in which the inflationary potential has a local minimum for a 
trapping of the field and a high temperature correction for a termination of infiation [21, 22]. 

In more recent studies of preheating processes [23-26] , primordial density perturbation 
spectra with blue tilt or, more generally, spectra having broad peak features at some fc* 
value are predicted in two-field models of infiation ending by preheating (going through the 



parametric resonance) [23] and even in one- field inflation models (in particular, in models 
with tachyonic preheating after small-field inflation [24, 25]). The characteristic values of A;* 
(i.e., the fe- values of modes which become dominant as a result of tachyonic amplification or 
parametric resonance) are different in different models. In preheating after chaotic infiation 
[16, 27] where the perturbations are amplified by parametric resonance, power spectra are 
peaked at scales k/aH ^ 1, whereas in models of small-field inflation the spectra may be 
peaked around the Hubble scale. In preheating after hybrid infiation, the typical scale of 
A;-values amplified by the tachyonic instability must be sub-Hubble if the phase transition is 
fast, i.e., if it takes less than a Hubble time [1, 26]. Clearly, the A;- value of the dominantly 
amplified mode is an important characteristic of a preheating model because large curvature 
and density perturbations of the Hubble size may result in an abundant production of PBHs 
and GWs. 

In the first stage of tachyonic preheating in models of hybrid inflation the ampliflcation of 
initial quantum fluctuations of the non-inflaton fleld is realized due to the classical inflaton 
rolling or, in a case of the small initial velocity of the inflaton fleld, due to processes of 
quantum diffusion [17]. The former case is characterized by the existence of a period of the 
linear evolution of the non-inflaton fleld. The evolution in this case can be studied using the 
cosmological perturbation theory [28-30] or 5N approach [31-34]. 

In the present paper we want to study in detail the dependence of the primordial curva- 
ture perturbations produced by hybrid inflation with tachyonic preheating on parameters of 
the inflationary potential. We show, in particular, that perturbation amplitudes strongly de- 
pend on the mass of the waterfall fleld m^. More exactly, it depends on the relation \m'^/H'^\, 
where H is a. value of the Hubble parameter during inflation. At small values of this relation, 
\m't/H'^\ ~ 1, the waterfall transition is rather slow, and the expansion of the Universe can 
not be ignored. Our main aim is to predict concrete values of the perturbation spectrum 
amplitudes, as well as a form of the /c-dependence of the spectrum. So, we preferred to use 
the results of the numerical (rather than analytical) solution of the equations for the time 
evolution of the waterfall field. The approximate analytic expression for the power spectrum 
amplitude, according to which one has, roughly, [8] 
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f ) , k<h, (1.1) 



was obtained in the limit of the fast waterfall transition, when |?7i?,| ^ H^, and can not 
be a priori used in a whole region of the parameter space. This formula works well when 
number of e-folds of cosmological expansion during waterfall is not too small (say, A^ > 0.1). 
If A^ ^ 0.1, eq. (1.1) must be multiplied on the factor which is proportional to A and can 
be much less than 1 [8]. 

The plan of the paper is as follows. In section 2 we calculate the time evolution of 
the power spectra of the waterfall field perturbations and the spectra at the end of the 
waterfall. The curvature perturbation spectra from the waterfall, as a function of model 
parameters, are calculated in section 3. In section 4 we give estimates for the possibility of 
PBH production in hybrid infiation model, with taking into account the non-Gaussianity of 
produced perturbations. Section 5 contains our conclusions. 



2 Calculation of waterfall field amplitudes 

We consider the hybrid inflation model which describes an evolution of the slowly rolling 
inflaton field (j) and the waterfall field X) with the potential [35, 36] 

V{ct>, X) = U^ - ^X^ + \m'4>' + \l<P\'- (2.1) 

The first term in eq. (2.1) is a potential for the waterfall field x with the false vacuum at 
X = and true vacuum at Xo = 2M^/vA = v"^ . The effective mass of the waterfall field in 
the false vacuum state is given by 
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rn^(<^)=7(<^'-0c), <fi^^^. (2.2) 



At (f) > (j)j. the false vacuum is stable, while at (j) < (pf. the effective mass-squared of x 
becomes negative, and there is a tachyonic instability leading to a rapid growth of x-modes 
and eventually to an end of the inflationary expansion. 
The evolution equations for the fields are given by 

4> + 3H,j)-V^^ = -^{m^ + jX^), (2.3) 

X + SHx- V\ = (2^/2M2 - ^d? - Ax')x- (2-4) 

Before the waterfall transition, i.e., at 0^ > </)^, the waterfall field is trapped at x = 0, so we 
can consider eq. (2.4) as an equation for the vacuum fiuctuation 5x- 

One should note that, as we will see below, the typical scale of the fluctuations amplified 
during the waterfall transition is of the order of the Hubble scale, or larger. In such a case, 
the back-reaction effects (due to perturbations of metric) may be non-negligible. The taking 
into account the back-reaction effects, i.e., an using of the perturbed FRW spacetime for the 
equation of motion for the scalar field (p at the stage when the second field, x, is absent, can 
be most simply done, as is well known (see, e.g., [37, 38]), in the spatially fiat gauge. In this 
case, the perturbed equation coincides with the unperturbed one in the limit (/> = 0. The 
smallness of (j) is provided by the slow-roll regime. So, to minimize the back-reaction effects, 
we will work in the vacuum-dominated regime, i.e., we assume that the vacuum energy M^ 
dominates, so that the Hubble parameter is effectively a constant. 






H"^ = Hl = ^^, (2.5) 



and one should consider only the case of small (j), satisfying the condition (j)'^ <^ MpH^. 

Another way to minimize back-reaction effects (and we will use it below) is to choose 
the gauge with uniform (/)-slicing, in which 6(j) = 0. 

As for eq. (2.4), the back-reaction effects are small because the unperturbed value of 
(x) is equal to zero [1]. 

The solution of eq. (2.3) (in which we ignore gradient term in accordance with our 
choice for the gauge) is (for t > tc, tc is the critical point when the tachyonic instability 
begins) 
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Figure 1. The numerically calculated dependence of the module of 5xk on time for two values of k: 
k = O.OliJc (solid curve) and k ~ He (dashed curve). Other parameters used for the calculation are: 
/3 = 100, r = 0.1, He = 10^^ GeV, (pc = O.lMp. For this set of parameters, the waterfall ends near 
Ht w 3.5, but we proceed the curves further just to show the asymptotic behavior. 



The scale factor a is normalized to 1 at the time of the beginning of the waterfall (at 
t = tc = 0), so 

a = e^^K (2.7) 

(2.8) 



The conformal time in this case is 



?y 



-Hct 



Hr 



aHr^ 



corresponding to the de Sitter expansion. 

From eq. (2.4), one obtains the equation for Fourier modes of 6x'- 

6xk + 3H6xk + (^ - f^H^c + 70') ^Xk = 0. 
Here, the parameter /3 is given by the relation 



(2.9) 



(2.10) 



Substituting the solution for (p in eq. (2.9) and introducing the new variable, u = a5x, 
one obtains the equation 



u'l + ik'^+f/ir]))uk=0, fJ.'^iv) 



^{\VH,\'' - 1) - 2 



(2.11) 



Here, primes denote the derivative with respect to conformal time ry. The normalization at 
early tim 
usual, is 



early times, when k ^ fi, is u = -r^e ''. The definition of the power spectrum of n, as 
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Figure 2. The numerically calculated spectra Vu{k) at different moments of time, for ji ~l,r = 0.1, 
He = 10" GeV, <?!)c = 3xlO-6A/p. Frombottomto top, -yyiJ = 10^, 10^, lO^^, lO^^, 10-^ lO"'^, lO'^. 
For these values of parameters, the waterfall ends at —rjH w 10^^. 

We show some results of numerical calculations in figures 1-3. In figure 1 we show how 
5xk depends on time for two different values oi k {k = H and k <^ H). It is seen that time 
evolution of both modes after the beginning of the waterfall is almost the same. Following 
[6] we assume that the growth era ends when the last term in right-hand side of eq. (2.3) 
becomes equal to the preceding one, i.e, when 



{{6xr) = ^ = xii- (2.13) 
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The asymptotic growth law of fe-modes of the waterfall field is approximately [3] 



Sxh-e^""', s = ^l + P-^. (2.14) 

It follows from eq. (2.14) that at large /3 the asymptotics of x is e^" , and the time scale 
of a period of the tachyonic instability is much less than a Hubble time. 

The typical example of the calculation of the power spectrum of u is shown in figure 
2. This Figure is similar with figure 3 of [3], where analogous quantities for the case /? = 
100, r = 0.1 were shown. It is seen from our figure 2 that even for the case /3 = 1 the waterfall 
is still effective, however it is much slower (it takes ~ 20 e- folds for /? = 1 while in the case 
of /3 = 100 the number of e-folds is ~ 3.5). 

Figure 3 shows the form of Vsxif^) at the moment of the end of the waterfall (determined 
by eq. (2.13)), for several parameter sets. Note also, that we impose an artificial cutoff of 
large-A: modes in our numerical calculation, which corresponds to considering only the water- 
fall field modes that already became classical at the beginning of the waterfall. Technically, 
the cutoff is imposed at the local minimum of 'P5i^(A;)-curve. 
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Figure 3. The numerically calculated spectra Vsxi^) ^^ the moment of the end of the waterfall. For 
all cases, r ^ 0.1, He = 10" GeV, 0c = O.lMp. From bottom to top, l3 = 1500, 100, 7. 

After a calculation of the power spectrum of 5x, we calculate the spectrum Prg^\2{k) at 
the end of the waterfall by the formula followed from the convolution theorem [39] : 



P, 



(sxy 



•ik) 



(2^)2 



d'k'Ps^{k')Ps^{\k-k'\). 



(2.15) 



3 Curvature perturbation spectrum 



The main equation for a calculation of the primordial curvature perturbation (on uniform 
density hypersurfaces) is [40] (see also [41, 42]) 

HSpnad 



c 



dt- 



P + P 



(3.1) 



where the non-adiabatic pressure perturbation is 6pnad = Sp — cj6p and the adiabatic sound 
speed is c^ = p/p. The formula (3.1) follows from the "separated universes" picture [31- 
34, 40] where, after smoothing over sufficiently large scales, the universe becomes similar to 
an unperturbed FRW cosmology. In our case, one has 



6p, 



'nad 



p 

SPx -- ^Px- 

p 



(3.2) 



Energy density p and pressure p is a sum of contributions of (p and x fields. Eq. (3.2) takes 
into account that in Spnad there is no contribution from (f) field. 

Everywhere below we will use for the fluctuations of the x field the more simple notation, 
than that is used in the previous section: x instead of 6x- It is just the same because (x) = 0. 

Our numerical calculations of Xk{t) in section 2 (two examples are presented in figure 
1) show that a form of the t-dependence of Xk weakly depends on k. Therefore, we assume 
that, approximately, the time dependence of x(x, i) is the same as of Xfc(i) at A; ~ fc*, where 
k-t is the peak value of V^^^. We designate this dependence (with arbitrary normalization) by 
f{t) (this function depends on the model parameters /3,r). Then, one has 



x(x,t) = C(x)/(t), 



(3.3) 



i.e., the dependence on a location is factorized, and 



x(x,t) = C(x)/(t) = x(x,t) 



m 
fit)- 



Expressions for the energy density and pressure of the x held are 

ml{t) 



Pxi^^t) 



^'"'x'lx,*) + -x'(x,t) - £2 ivx(x,t)r 



mlit) 



Pxi^^t) = ^^xH^,t) + ^X'(x,t) + ^ |Vx(x,i)'' 



1 

1 
2^2 



(3.4) 

(3.5) 
(3.6) 



(where, for brevity, m^ {4'i't)) = ''T^'xi't))- 

For calculations of the curvature perturbation spectra we will use the spatially averaged 
energy density and pressure, i.e., we do in equations (3.5, 3.6) the substitutions x^ ~^ (x^); 

x'^(x'>, IVxP^(IVxP). 

We proved the relative smallness of gradient terms in these expressions using the ap- 
proximate method suggested in [6] for the case when cosmological expansion during the 
waterfall is ignored. Namely, due to a peak feature in P^(A:)-dependence (see figure 3) one 
has, approximately. 



(ivx(x,t)r 



d-'kk'\xk\' ~ K{t){x'{t)). 



(3.7) 



For values of parameters we operate with. A;* ~ He- At the end of the waterfall {t = tend), 
the estimate for the gradient term will be 



(ivxr 



u2 2 



(3. 



and the relation between gradient and time derivative terms in (3.6) (we denote it with letter 
D) can be written as 



D 



(ivxr 



1 



a2(x2) ' sMtendV' ^^'^^ 

where we have used the asymptotic behavior (2.14) for this estimate. For P ^ 1, D is 
negligible due to the large factor s, s » 1. For /3 ~ 1, s ~ 1, but a » 1. We conclude 
that D is always small and the contribution of spatial gradient in (3.5, 3.6) can be always 
neglected (note that in case when /3 is small, the peak in 7^5^(A:), as seen from figures 2 and 3, 
is very broad, so in this case the smallness of gradient terms is proved by the straightforward 
calculation of (3.7)). 

For a calculation of 5pnad we need expressions for fluctuations of energy density and 
pressure. They are given by 
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5x\ 



(3.10) 



(3.11) 



where Sx^ = X^ ~ ix^)- Averaged values of x^ and x^ are calculated taking into account the 
equation (2.13): 



ix'it)) = xlr ^'^'^ 



{r{t))=x 






(3.12) 
(3.13) 



Here, tend is a moment of time of an end of the waterfall. Finally, for the average values of 
the energy density and pressure, one has 
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The values of p^ and p^h are 

P<l> -- 



vi<p) + -0^ 



(^^ ] x\ 

Uitend) ' ^"' 



vi^) + \^\ 



(3.14) 

(3.15) 
(3.16) 



where, for our case, V{(l)) = g^^i^c/)^. 

For the curvature perturbation, we have the integral 



c = c. 



Hrdt 
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The relation between curvature perturbation and x^-spectra can be written as 

V( = A Vsx^itend), 



(3.17) 



(3.18) 



where the value of the spectrum in the right-hand side is calculated at a time of the end of 
the waterfall, tend, and A is to be determined. 

It is easy to check that for a calculation of the spectrum 'P^^a in equation (3.18) one 
can use the expression for Vy2 given by the convolution formula (2.15). 

The time dependence of x^(t) needed for the calculation of the integral in (3.17) is given 

by 

X'(t) = X'(tend) (y£^) . (3.19) 



Using this relation, we can extract the value of A: 
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(3.20) 

We show results of the calculation of V^ik) in figure 4. It is seen that the spectrum can 
reach rather large values (of order of 1) for the case of /3 ~ 1 (if r is in a broad interval, say, 
0.1^0.001). 
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Figure 4. The numerically calculated spectra 7^^(fc) at the moment of the end of the waterfall. 
Parameters used for the calculation are: r = 0.1, He = 10^^ GeV (all curves); (f>c = 3 x 10~^Mp (for 
13=1,2, 7); (pc = O.lMp (for (3 = 100, 1500). 



In figure 5 we show the model parameter regions that lead to rather large values of 
Vt^. We also show, on the same Figure, the parameter region where waterfall is not effective 
(classical regime is not reached) and boundary of (/>c = 1 (in units of Planck mass Mp). 

At the end of this section, we numerate the constraints on the parameter space which 
follow from the assumptions used for an obtaining of the main results. 



1. Slow roll of the inflaton field: 



m^ < Hi 



(3.21) 



2. False vacuum dominance: 



M^ > -m^(f)i. 



(3.22) 



3. Small (/>-condition (for justifying of the omission of high powers of (f) in the inflaton 
potential): 

(t)c < Mp. (3.23) 



4. Small effects from back-reaction (slow-roll condition): 

(/)2 < MJ>H'^. 



(3.24) 



5. Condition for 6xk to be classical field [1] 



^""W 



(3.25) 



6. Condition for a dominance of the scenario with the inflaton's classical rolling [26]: 



<j)c > -yXv"^, or 7^/^ < 



(3.26) 
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Figure 5. The allowed ranges for the values of 7 and A and the corresponding maximum values of 
■p^(fc) that are reached. Left panel: r = 0.1 (to k, 0.5Hc). Right panel: r = 0.001 (to w 0.05i?c)- For 
both panels, ifc = 10" GeV. 

4 PBH production from non-Gaussian perturbations 

A production of PBHs (about these objects, see, e.g., reviews [43, 44]) during reheating 
process had been considered in works [20, 23, 24, 45, 46]. The classical PBH formation 
criterion in the radiation-dominated epoch is [47] 

6>6c- 1/3, (4.1) 

where 5 is the smoothed density contrast at horizon crossing. The Fourier component of the 
comoving density perturbation 5 is related to the Fourier component of the Bardeen potential 

^ as 

2 / k ^^ 



For modes in a super-horizon regime, ^^ ~ —(2/3)TZk ~ — (2/3)Cfci so (4.1) can be translated 
to a limiting value of the curvature perturbation [48], which is 

Cc ~ 0.7. (4.3) 

It is seen from (3.17) that the curvature perturbation generated by waterwall field has 
a negative sign, so, naively, the threshold (4.3) can't be reached. However, the perturbations 
must be considered with respect to the average value, so 

C^Co = C-(C), (4.4) 

and the condition for individual PBH formation is Co ^ Cc while PBHs will exceed the 
currently available limits on their average abundance already for | ((") | only slightly exceeding 
Cc (i-e., |(C)| — Cc ^ 1) [8], so, practically, the PBH constraint on ( is just |(C)| < Cc- 
For the distribution of Co we may write 

Co = -{9' -{9')), (4.5) 



10 



where g is gaussian (in our case, g ~ 5x), and for the average of C,^, using known properties 
of Gaussian distributions, 

{Cl) = 2{gY = 2m\\ (4.6) 

so, in terms of (^q) the significant PBH formation will happen if 

(Co') = /n(^)f > 2d-l- (4.7) 

In case of P^-spectra shown in figure 4, the curve corresponding to /3 = 2 is close to 
the bound (4.7) while the curve for /? = 1 has (Cg) ~ 7 which makes that set of parameters 
forbidden by the PBH formation constraint. 

The mass of the PBHs produced from curvature perturbation spectra presented in figure 

4 can be estimated as follows. The horizon mass at the end of inflation is 

M, .. {H~^fp = ^ ~ 10^ g, (4.8) 

and using a well-known dependence for the horizon mass corresponding to fluctuation having 
wave number k, Mh ~ k~'^ (see, e.g., [49]), we have 

Mbh ^M,^M, ^^^hl^^ ' . (4.9) 

For the case of figure 4, curve for /3 = 1 corresponds to Mbh ~ 10^^ g (the mass range 
of non-evaporating PBHs that can constitute dark matter) while /3 = 2 corresponds to 
Mbh ~ 10^'^ g (such PBHs have already evaporated, but the products of their evaporation 
are, in principle, observable). 

For other values of inflation energy scale, or reheating temperature, which is connected 
to He by the relation [49] 

/90M2ff2N'/4 
Trh = 2^^ , g.- 100, (4.10) 

we sketch the expected range of characteristic masses of PBHs, that can be produced (see 
figure 6). The shaded region in the figure corresponds to large values of /? parameter (1 < 
/3 < 2). The results shown in this figure must be considered together with the particular 
cosmological constraints on PBH abundance. We leave this question for further studies. 

5 Conclusions 

We carried out numerical calculations of a contribution of the waterfall field to the primor- 
dial curvature perturbation (on uniform density hypersurfaces) C, which is produced during 
waterfall transition in hybrid inflation scenario. The calculation is performed for a broad 
interval of values of the model parameters. 

We did not consider the contribution to ( from inflaton rolling, which is dominant at 
cosmological scales. One should note that the simple quadratic inflationary potential used 
in the present paper can be easily corrected (see, e.g., [50]) to give a red-tilted spectrum at 
cosmological scales (converting, e.g., the original hybrid inflation model to a hilltop model 
[51, 52]), without any modiflcation of our predictions for small scales. 
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Figure 6. The range of characteristic masses of PBHs tliat can be produced by a hybrid inflation 
waterfall, as a function of inflation energy scale or reheating temperature. 

Main results of the paper are shown in figures 4, 5, 6. One can see from the figure 4 
that peak amphtudes of P^ strongly depend on the value of /3 (curiously, 7^^ ~ 1 corresponds 
to /^ ~ 1 in a broad interval of 7 and r). The peak values, k^, for small (3 are far beyond 
horizon, so, the smoothing over the horizon size will not decrease the peak values of the 
smoothed spectrum. Furthermore, the spectrum near peak remains strongly non-Gaussian 
after the smoothing. Our calculations based on the quadratic inflaton potential show that 
for (3 < 100 and in the broad interval of r the peak value k^ can be estimated by the simple 
relation: 






~ e 



-N 



(5.1) 



where A^ is the number of e-folds during the waterfall transition. The similar estimate is 
contained in the recent work [8] . 

We conclude that the strong growth of amplitudes of the curvature perturbation spec- 
trum at /? — )• 1, which had been anticipated in pioneering works [53, 54], really takes place. 
However, the condition which is very often used as a constraint, that the bare mass-squared 



of the x-field. 



-mx 



2\/AM^, must be much larger than H^ (the so-called "waterfall con- 



dition" [55]) seems to be too restrictive. Our results, presented in figures 4 and 5, show that 
only at — w,?, ~ H^ (when /3 ~ 1) the amplitude of the curvature spectrum Vc_ becomes close 
to one, i.e., enters a region that can be constrained by PBH data. In figure 6 we show the 
region of PBH masses which can be produced due to the waterfall transition in hybrid infla- 
tion, in the case of large (3 parameter. Abundances of these PBHs should be constrained by 
data of PBH searches. It would be very interesting to analyze the corresponding constraints 
following from data of relict GW searches (see, e.g., [26, 56, 57]). 
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